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Gauge invariant quintessence perturbations in the quintessence and cold dark matter (QCDM) 
model are investigated. For three cases of constant equation-of-state (EOS) parameter, linear scalar 
field potential, and supergravity scalar field potential, their perturbation evolutions have a similar 
dependence on EOS parameter and scale, but they have different sensitivity to the initial conditions 
due to the different shapes of the quintessence potential. They have a minor effect on primary CMB 
anisotropies, but change the secondary CMB effect in different ways. The first case is insensitive to 
initial quintessence perturbations and only modifies the Integrated Sachs- Wolfe (ISW) effect within 
a factor of 2. The other two cases are sensitive to initial conditions at large scales and could affect 

■ the secondary CMB anisotropies drastically depending on how smooth the initial perturbations 

are. This makes it possible for future cosmological probes to provide constraints on quintessence 

, properties. 

(N 

^ , I. INTRODUCTION 

i Recent observations of an accelerating universe [J-Q-Q imply the existence of dark energy characterized by a negative 
fNJ ' pressure to density ratio, also known as the EOS parameter w. There have been discussions of the cosmological 
constant, the dark energy evolving according to a specific EOS parameter 0Q|i an d the dark energy consisting of 
a dynamical cosmic scalar field, the quintessence lllj. In contrast to the cosmological constant, the quintessence 
component in QCDM model has kinematic behavior and can develop perturbations. It has important features of time- 
evoluting EOS parameter and scale-dependent effective sound speed [T3. IlH ] . One significant issue is the sensitivity 
to initial perturbation conditions, which attracts much attention. 

For some QCDM models, the quintessence potentials can be approximated with constant EOS parameter, and the 
perturbations evolutions are insensitive to the initial quintessence perturbations, with minor effect on observations 
UH [Hi Hdll < For exponential potential V(Q) = Ve~^ c ^ M ^ H,0,E3j although quintessence perturbations may stay 
non-zero at some time for different background attractor solutions, they die out in matter-dominant time for physically 
reasonable models which shows the insensitivity again. For the tracking quintessence models (2£j , large-scale 
Q H i non-adiabatic perturbations can grow or stay constant before entering tracker regime, but have to get suppressed 
afterwards |2j|. However, the first-order matrix formulation developed latter [22j applies to general quintessence 
potential cases on all linear scales and indicates the possibility of non- vanishing entropy perturbations. Recently some 
other QCDM models with special quintessence potentials, such as linear V(Q) [2(j and supergravity V{Q) [25|. lead 
C$ • to special background evolutions. So it will be interesting to see the perturbations behavior and their sensitivity to 
initial conditions in these models, as well as their effects on cosmic microwave background (CMB) jj| and baryon 
• i— i , mass power spectrum [j| in the universe. 

In this paper, we describe the universe by a simple QCDM model with the conformal-Newtonian metric in section 
2. The cosmological perturbations are solved for constant w cases in section 3, for \me&r-V(Q) cases in section 4, 
and for supergravity-y(Q) cases in section 5. The results are analyzed by the matrix formulation and their effects on 
cosmological observations are discussed in section 6. 
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II. QUINTESSENCE AND COLD DARK MATTER (QCDM) UNIVERSE 

Here we consider the uniform matter and quintessence field background with small perturbations in a perturbed 
flat FRW metric. Focusing on scalar perturbations, we can use the flat conformal-Newtonian metric for the gauge 
invariant approach [2^.l2^ |: 

ds 2 = a 2 (r]) [-(I + 2<f>(r),x))dr] 2 + (I - 2<f>(r),x))5 ij dx i dx j ] (I) 

where rj is the conformal time and \$>(rj, x)\ <C I is the gravitational potential. We define ' = J--, H = 2- as the "Hubble 

parameter" in the conformal-Newtonian version, 770 as today's conformal time, and a = as the normalized 

scale factor. For the quintessence scalar field, we assume the standard Lagrangian Cq = —hQ^Q'^ — V(Q) and 
decompose the quintessence field as uniform background with small spatial perturbations: Q(j), x) = Qo(ri) + 5Q(r], x), 
\SQ\ <C \Qo\, (5Q(i], x)) — 0, where ( ) means global spatial average. We assume that the quintessence part and matter 
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part only interact with each other gravitationally. So the wave equation of the quintessence field — + dV/dQ = 
gives its background and perturbative components: 

QZ + 2HQ' Q + a 2 V Q (Q ) = 0, (2) 
SQ" + 2HSQ' - V 2 5Q + a 2 V QQ 6Q - 4Q' 1>' + 2a 2 V Q ^> = . (3) 

The energy-momentum tensor of this quintessence field is Tq U = Q 4i Q M + {—\Q, a Q ,a — V) g^ v ' . 

For the matter part, we assume its energy-momentum tensor has the form of the perfect fluid: T Miy = (p m + 
p m )U^U v +p m g» v , where = dx^/dr = [(1 - $) /a, V l /a] is the 4-velocity of the fluid, V 1 = dx^drj = a dx^dt < 1 
is its 3-velocity, p m (rj,x) = (p m ){rf)[l + 6 m (r),x)], and \5 m (rj,x)\ <C 1. If we neglect all pressure effect (\p m (j],x)\ -C 
Pm(Vjx)i |V<5p m | <C / o m |V<£>|) on scales which are not too small (> 10 Mpc), the energy-momentum-conservation law 
gives background and perturbation components after some manipulation: 

(p m ) = , (4) 

C + H6' m - 3$" - 3ft $' - V 2 $ = . (5) 

The matter perturbation 5 m can be formally obtained in terms of a and <£>(£;, ?y) in fc-space, where k is the comoving 
wave- vector (— V| = k 2 ). 

To complete this system, we consider the Oth-order 0-0 component and the lst-order i—j component of the Einstein 
equations: 



m 2 

— 7T- = 8wG 



d 3 2a 



, - V(Q ) 



(6) 

$" + 3ft$' + (2W' + H 2 )<S> = AnGa 2 ( - + - V^q) . (7) 



Given the initial conditions, V(Q), and other parameters, eq.J3ll|) give the background Qo(r]) and a(rj), and eq.J31 
0) give the perturbations <5Q and Due to the convenient conformal-Newtonian metric, the matter perturbation 8 m 
drops out and can be found later via eq.©. 



III. CONSTANT w CASE 



Here we follow the EOS parameterization in ^| to study the general properties of the quintessence field instead of 
some specific quintessence field potential V(Q). 



w{rj) 



(pg) _ % v (Qq) 

( PQ ) Q£ + V (Q ) 



Then we have 



2cr 1 — w 



1 + w 



(pq)- 



By manipulating eq.(|21IHJ), we can express Vq and Vqq in terms of the EOS parameter w(rf): 
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With eq. (JTO| EJ and definitions 



Qo = a iVo)^J (1 + w){pq)(t} Q ) i//, ^'(r ?0 ) = l, 
<5Q ee a(r] )^ (1 + w)(p Q )(rj ) Sip, 



(8) 

(9) 

(10) 
(11) 

(12) 
(13) 
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we can put the quintessence wave equations eq. EI) in a more convenient form: 



Si//' + (2H + )Si// + <-V 2 -^-(l-w) 
1 + w ! 2 



a „ i9 . 7 3w , 

— -w 2 ~ + ^- 

a 2 2 



3H(l-iy) 



1 + w 



V>'$ = 0. 



3Wu/ ^ 



(14) 



(15) 



In this formulation, any function V(Q) is equivalently described by the corresponding EOS parameter w(rj), which 
can be much more easily specified and compared with other dark energy models. One simple example would be 
w = constant, which has the exponential- like potential V(Q) shown in Figure ^ 
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FIG. 1: Scalar field potential V(Q) for constant-u; cases. The curved line is for w = —2/3. The scalar field rolls from top-left 
corner down to the today's value indicated by the dot. The horizontal line is for w = — f (AC DM). 



In this special case, ea . (|9" HT2l II 41) give the solution for the quintessence field background: 



(pq) 



(pq(vo)) 



fl3(l+tu) ■ 

Then the previous closed system of a, Q' , $ and 5Q can be reduced to a system of a, $ and Sijj in /c-space: 



T-i 2 

% = (Hoaof 



n, 



a 3 a 3 ( 1 +™') 



$" + 3ft$' + (2W + H 2 )$ = -(iJ ao) 2 (l + w)^Q 



<ty" + TKH' + \k 2 - -(1 - w) 



a 2 2 



q1+3«j g(l+3tu)/2 2 y a( 1+3tu )/ 2 
^ - 4 , - 37^(1 - w) ,_, . „ w » 

q(1+3uj)2 v ' q(1+3u))/2 



0, 



(16) 

(17) 

(18) 
(19) 
(20) 



where Q m and £Iq are the density parameters defined as usual. If w = — 1, the quintessence perturbations vanish in 
ca. (|19fl and only its background contributes. This special case thus reduces to the AC DM model. At early times, this 
model can be well approximated by the flat CDM model (Einstein-DcSitter space). Then ea. (|18j) gives a(ry) oc rf , 
U ~ 2 - 1 and eq.JEJ gives <I> independent of time for the growing modes of matter perturbations S m oc a(rf), as usual. 
At later times, the quintessence part begins to become important and the solution of a(rj) from eg. 1181) will deviate 
from n 2 more and more as time goes on. So the coefficient 2Ti' + 7i 2 in eq. H19|) is not zero any more and makes <& 
decay. This is what causes the ISW effect in low-Z multipoles of the CMB power spectrum. 

Recent cosmological observation of the accelerating universe requires w(r]o) < —1/3. The w < — 1 cases would 
correspond to negative kinetic term in the standard quintessence Lagrangian or other more complicated models, 
which we won't discuss here. So from now on we just consider —1 < w < —1/3 in all our QCDM models. The initial 
conditions are a(0) = 0, |$(0)| ~ 10~ 5 and $'(0) = for modes outside horizon, based on inflation theory [33]. For 
the same reason, |<5-0(O)| ~ |$(0)|. Because in ea. ipU)) the restoring term (Vqq oc H 2 ) is bigger than the damping 
term (oc H) which is bigger than the driving term (oc Hd~^ 1+3w ^ 2 ) at early times, Sif> with any initial value which is 
not too big (|<5-0(O)/$(O)| < 10 4 ) will oscillate and get damped quickly without affecting the $ evolution before the 
last-scattering epoch. So the primary effects on the CMB won't get modified much and are insensitive to the initial 
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condition on the quintessence perturbations. This is one important result within many QCDM models llallfj- Here 
we can see that it is mainly due to the huge Vqq, i.e., the huge quintessence mass at early times. For now we just 
take the smooth initial condition 5ip(0) = 5ip'(0) = for convenience. Let us look at the behavior of this system at 
various spatial scales. 

At small scales (k ^S> 1), by ea. (|20|l the growth of Sip is strongly suppressed and stays negligible compared with $. 
The effective sound speed of the quintessence perturbations approaches 1 (c 2 = Spq/5pq ~ 1). Then we can drop the 
5ip terms in the above system and the evolution of $ only comes from background effect without a dependence on k 
in this region: 

+ m& + (2W + H 2 )<$> = -^(H a ) 2 (l + ■ 

This system behaves like the ACDM model with a different EOS parameter w and the additional r.h.s term. If uu is 
bigger the dark energy takes over dominance earlier and a deviates from n 2 earlier too, which makes <£> decay more. 
In addition, the negative r.h.s brings down <E> more. So compared with the ACDM model, small-scale quintessence 
perturbations in the QCDM model do not grow and only the background effect contributes, which drags down the 
gravitational potential more and corresponds to a bigger ISW effect. 

At scales comparable to or larger than horizon size (k < I), the quintessence perturbations dip can grow to the 
order of the gravitational potential $ by ea. H2Q(l with smaller effective sound speed (c 2 < I). Then it will back-react 
on the <f> evolution via ea. (|I9|l . Because of the positive coefficient, the quintessence perturbations tend to raise the 
gravitationally potential $ to make it decay less. This is different from the k ^ 1 case, and agrees with the result 
of a phenomenological approach |23||. Thus in this QCDM model the evolution of perturbations and gravitational 
potential depends on scale. 

With definition on horizon scale ao = numerical calculation gives $>(r]) and 5ip(n) once we specify the parameters 
k and w. Corresponding to the assumptions we made before for this linearized system, we only consider k < 800. For 
w = —0.5 and k = {0.01, 5, 500}, the time evolution is shown compared with the w = — 1 case in Figure^ 
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FIG. 2: Time evolution of the gravitational potential &(n) for the constant-w case. The three thin lines correspond to 
k — {0.01, 5, 500} from top to bottom for the w — —0.5 case. The thick line is the scale-independent w = — f case (ACDM). 



IV. LINEAR V(Q) CASE 

Linear scalar field potential was looked by Andrei |32( i n inflation theory long ago. Recently it is reviewed by 
Dimopoulos and Thomas |26| as a quintessence field within a quantum mechanically stable action. Due to the large 
Z-moduli of the wave function factor, the higher order potential derivatives d n V/dQ n (n > 1) can be neglected and 
only the linear term remains. After shifting the origin and considering the symmetry of Q <-> —Q, the typical V(Q) 
can be written as: 

V(Q) = -Vi Q, Vi = constant > 0. (21) 

This linear potential may eventually lead to negative total energy density and result in a rapid crunch. Let us consider 
the closed system defined in §2 to study the behavior of both background and perturbations for this model in detail. In 
order of magnitude, Eq.@ gives Q' ~ a^oVi ~ a Vi/H and eq.@ gives Q' 2 /2a 2 ~ \ViQ\ ~ (pq) ~ 3Hq£Iq/8itG for 
1 + w ~ 1. Combining these two, we will have Vi ~ Hq^^q/SttG ~ H 2 m p (h = c = 1) and Q ~ ^3Qq/8ttG ~ m p . 
So this simple analysis sets the range of magnitude for the quintessence field quantities. Now we can redefine the 
quintessence field for later convenience: 



1 1 = H " X '^G ^' Q = V IS§ ^' ^ ~ ^ ~ ^ = constant > °- ( 22 ) 
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The closed system of equations can then be rewritten as: 



Qo + 2HQ' - H%ala 2 V x = , 



ttt — H ao 



\ 



Q'o 2 



-ViQc 



SQ" + 2HSQ' + k 2 SQ - AQ' & - 2H$a 2 ) a 2 V 1 <f> = , 



$" + 3H$' + (2H' + H 2 )<S> = -n Q 



(23) 

(24) 

(25) 
(26) 



where the perturbation equations are in fc-space. For the background, we take the initial condition as a(0) = 
and Qb(0) = based on inflation theory. The Qo(0) value is chosen by requiring that today's quintessence energy 
density fits the observation, pq(tj ) = 3Hq$}q/8itG, or equivalently, Q' 2 j2H 2 d^ — ViQq(t]o) = 1. Then this describes 
an initially static field rolling down a potential slope. Its effective EOS parameter evolves from -1 to some value 
— 1 < w(r)o) < today. Different slope values V% correspond to different Qo(0) and EOS parameter w(r]), but the 
same energy densit y pp (?7n) today. The correspondence is listed in Table ||| The detail has been worked out in the 
"doomsday" model j22|. F° r the w(r]o) = —0.5 case, the potential V(Q) is shown in FigureOl 



TABLE I: Quintessence potential slope and starting position corresponding to w(r]o) 

Vi Qo(0) w(r]o) 
2.738 -0.635 -0.5 
2.424 -0.657 -0.6 
2.077 -0.696 -0.7 
1.677 -0.774 -0.8 
1.172 -0.980 -0.9 




FIG. 3: Linear scalar field potential V(Q) for the w(j]o) = —0.5 case shown as the slope. The two dots are the starting and 
today's positions of the field variable from left to right. The horizontal line is the w = — 1 case (AC DAI). 

For the perturbations, again we take the initial condition as 5Q'(0) = $'(0) = and |$(0)| ~ 10~ 5 , based on 
inflation theory. Ea. 1)25(1 describes the oscillation of an intitially static quintessence perturbation field driven by the 
gravitational potential. At small scales (k ^> 1) the quintessence perturbations stay small. So the $ evolution only 
depends on the background. The behavior of the perturbations is not sensitive to k or 6Q(0), as in the constant-w 
case. However, at scales comparable or bigger than the horizon size (k < 1), due to the zero quintessence mass 
(Vqq = 0) in the restoring term, the quintessence perturbation will stay near its initial value with a shift driven 
by the gravitational potential. Then the whole system becomes sensitive to the initial condition of quintessence 
perturbations, which is quite different from the constant-w case. For the w(r]o) = —0.5 special case, the dependence 
on k and 5Q(0) is indicated in Table ILTI 

If we take the smooth initial condition SQ(0) — 0, then in terms of different k the perturbation evolution result 
is shown in Figure Again we see that the scale dependence of perturbations evolutions is similar to the constant- 
w case in the last section. The quintessence perturbations stay suppressed, corresponding to more decay of the 
gravitational potential at small scales, but grow to the order of the gravitational potential accounting for less decay of 
the gravitational potential, with the transition scale around k = 5. Besides, the decay of the gravitational potential 
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TABLE II: Dependence on initial perturbations in the linear-V(Q) case, for w(r]o) = —0.5 





k = 


0.1 


k = 


1 


k = 


10 


dQ(0)/<P(0) 


dQ(ri )/<P(0) 






$(t7o)/<E>(0) 


oQ(-no) mO) 


$(r?o)/$(0) 


10 


12.077 


3.862 


9.411 


3.262 


0.147 


0.741 


1 


1.752 


1.116 


1.469 


1.054 


0.108 


0.730 





0.605 


0.811 


0.586 


0.808 


0.104 


0.729 


-1 


-0.542 


0.506 


-0.296 


0.563 


0.099 


0.728 


-10 


-10.867 


-2.240 


-8.240 


-1.646 


0.060 


0.717 



is smaller than the constant-w case because its quintessence EOS parameter stays close to -1 at early times and only 
rises to w(r)o) recently, so that the dark energy emerges later. 

l 

0.9 

o 0.8 

o 

o 0-7 
. 6 
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FIG. 4: Time evolution of the gravitational potential for the \meax-V(Q) case with w(rio) = —0.5 and smooth initial 

condition. The three curves correspond to k = {0.01, 5, 500} from top to bottom. 




V. SUPERGRAVITY V(Q) CASE 

The scalar field potential can also come from extended gauged supergrayity, which has a close relation to M/string 
theory and extra dimensions. Masses of ultra- light scalars in these models (2£| are quantized by the Hubble constant: 
to 2 = n 2 H 2 . If the de-Sitter solution corresponds to a minimum of the effective potential, the universe eventually 
becomes de-Sitter space. If the de-Sitter solution corresponds to a maximum or a saddle point, which is the case 
in all known models based on N — 8 supergravity, the flat universe eventually stops accelerating and collapses 
to a singularity. For example, all known potentials of the gauged N — 8 supergravity have the universal feature 
n 2 = {—6,4, 12}. Along the tachyon direction n 2 — —6, the scalar potential can be written as: 

07TG 



Vq = constant > 0, Q = — , 



rn p V 87rG 

where Vq and Q are dimensionless and in the order of 1. Again we apply the closed system defined in §2 to this model. 
The background and perturbations equations are: 

Qg + 2HQ' - 6H 2 a 2 n Q V a 2 Q = , (28) 
a' 



¥ ~ V i 5 6s 2 " 



H 2 a 2 Q Q V Q Qo(l - Ql) , (29) 



SQ" + 2HSQ' + {k 2 - 6H 2 a 2 a 2 n Q V )SQ - 4Q' & - 12H 2 a 2 a 2 n Q V Qo<S> = , (30) 
$" + 3H$' + (2H' + Tt 2 )$ = -i(?o 2 $ + \q' 8Q' + 3H 2 a 2 a 2 n Q V a Q a 5Q . (31) 

For the background, we again take the initial condition as a(0) = and Qd(0) = 0. We choose Qq(Q) = 
{0, 0.2, 0.3, 0.32} and adjust Vq so that today's quintessence energy density fits the observation, pq(t]o) — SH'qQq /8ttG, 
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or equivalently, Q' 2 / 'QH^o^VLq + Vq[1 — Qq(vo)] = 1- Then this describes an initially static field rolling down a curved 
potential slope. Its effective EOS parameter evolves from -1 to some value —1 < w(rio) < today. Different initial 
positions Qo(0) correspond to different Vq and EOS parameter w(r]), but the same energy density pq{t]o) today. The 
correspondence is listed in Table ITTT1 For the Qq(0) = 0.3 case, the potential V(Q) is shown in Figure [S] 

TABLE III: Quintessence initial position corresponding to Vo and w(t]q) 



Qo(0) 


Vo 


w (770) 





1 


-1 


0.2 


1.124 


-0.92 


0.3 


1.441 


-0.67 


0.32 


1.645 


-0.47 




FIG. 5: Scalar field potential V(Q) from supergravity for Qo(0) = 0.3 shown as the curve. The two dots are the starting and 
today's positions of the field variable, from left to right. The horizontal line is the w = —1 case (ACDM). 

The perturbations are analyzed as before. We take the initial condition as 5Q'(0) = $'(0) = and |$(0)| ~ 10 -5 . 
Eq.(j20J describes the oscillation of an initially static quintessence perturbation field driven by the gravitational 
potential. At small scales (k ^> 1) any initial quintessence perturbations get damped to be small, and the $ evolution 
only depends on the background. This system is not sensitive to SQ(0), as in the linear- V(Q) case. At scales 
comparable to or bigger than the horizon size (k < 1), due to the tachyon quintessence mass (Vqq < 0) in the 
restoring term, a quintessence perturbation will stay near its initial value and even grow, with a shift driven by 
the gravitational potential. Then this system is more sensitive to the initial quintessence perturbations than the 
linear-V(Q) case. For the Qo(0) = 0.3 case, the dependence on k and 5Q(0) is indicated in Table llVl 



TABLE IV: Dependence on initial perturbations in Supergravity- V(Q) case, for Qo(0) = 0.3 





k = 


0.1 


k = 


1 


k = 


10 


5Q(0)/$(0) 


5<3M/$(o) 


*(ijd)/*(0) 




<&M/<&(o) 




*fa))/*(0) 


10 


23.412 


3.016 


18.640 


2.589 


0.207 


0.753 


1 


3.018 


1.011 


2.521 


0.967 


0.159 


0.746 





0.752 


0.788 


0.731 


0.787 


0.153 


0.746 


-1 


-1.514 


0.565 


-1.060 


0.606 


0.148 


0.745 


-10 


-21.908 


-1.440 


-17.179 


-1.016 


0.100 


0.739 



This potential is unbounded from below and the theory is unstable. Since the potential can remain positive for 
|<3?| < 1 for a time longer than the present age of the universe before it finally collapses, this model is sufficient to 
describe the background evolution of the universe. However, the tachyon mass makes the quintessence instability 
develop and affect the gravitational potential, and thus the system becomes sensitive to the initial conditions at large 
scales (k < 10). This puts more constraints on the initial quintessence conditions and this supergravity model itself. 
For the smooth initial condition 5Q(0) — 0, this model has a similar scale-dependent perturbation evolution as the 
lmear-V(Q) model in last section. For Qq(0) — 0.3 and k = {0.01,5,500}, the time evolution $(77) is shown in 
Figure 
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FIG. 6: Time evolution of the gravitational potential for the supergravity-V(Q) case, with Qo(0) = 0.3 and smooth initial 
condition. The three curves correspond to k — {0.01, 5, 500} from top to bottom. 



VI. DISCUSSION 



In another viewpoint, we can apply the first-order matrix formulation in |22j to analyze the different quintessence 
perturbation behaviors in above QCDM models. At large scales (k <C 1), the relative entropy perturbation S and 
intrinsic entropy perturbation T |l7l Il8l Il9| have coupled evolution equation 



9 [ S \ _ „ / w Q ~Wf + ifSlf(wf - <? sQ )h 7/^/(1 





d\n(a/a ) \ V J \ -7/2 wq 



where / means the perfect fluid component, 7^ = 1 + Wi, and c 2 q = p/p. For most of the time in our QCDM 
models, CDM dominates and have Wf ~ 0, 0/ ~ 7/ ~ 7 ~ 1, and —1 < wq < —1/3. Then the eigenvalues are 

n± = I 2wq — c 2 s q — i ± J (c 2 s q + |) 2 — 2 . For the constant w case, wq — c 2 s q leads to Re(n±) = |(wq — 1/2) < 0, 
corresponding to decaying entropy perturbations. While for the linear V(Q) case, c 2 q ~ —2 leads to n± = (3wq + 
3,3u>q + |), corresponding to at least one growing entropy perturbation mode. This fits our previous analysis well. 

In all QCDM models that we studied above, the quintessence perturbations have the common characteristic feature 
of scale dependence. At small scales the perturbations do not grow and the gravitational potential decays more only 
due to background effects; at large scales the perturbations grow and sustain the gravitational potential. On the 
other hand, because of the different quintessence potentials, the three QCDM models have different sensitivity to 
the initial conditions on quintessence perturbations. The constant-u; model has a large quintessence mass at early 
times and damps perturbations fast so that it is insensitive to initial conditions at all scales. In this model the 
quintessence has no effect before the last-scattering time. In contrast, the linear-y(Q) model and supergravity-V r (Q) 
model have zero and tachyon quintessence mass respectively, and their perturbations remain near their initial value 
(<5<3(0)/<I>(0) < 10 4 ) at large scales. Since the EOS parameter w in these two models approaches -1 at early times, 
their dark energy dominates only very recently, like the ACDM model. So the total back-reaction of dark energy 
perturbations (background times relative perturbations) on other perturbations remains negligible before the last- 
scattering time, when the universe evolves like a oc rj 2 , Q ~ Qo(0), Q' — 0, SQ ~ SQ(0), and $ ~ $(0). So they do 
not change the primary CMB anisotropics either. But later when the dark energy background become dominant, this 
back-reaction of dark energy perturbations will affect other perturbations drastically, depending on how smooth the 
initial perturbations are. 

For smooth initial conditions in all models the dark energy has no primary CMB effect except that its background 
contributes to the distance to the last-scattering surface. But later on both the dark energy background and per- 
turbations will contribute to the change of the gravitational potential, which may lead to a considerable secondary 
CMB effect such as the ISW effect [3(|. For the scale-invariant primordial spectrum P$ = |$(0, k)\ 2 = A 2 /k 3 , the 
large-scale CMB power spectrum due to the SW effect is: 



C sw = ^ 2 K 2 
1 36tt 2 Z(Z + 1) 1 ' 



f° 

Kf - 21(1 + 1) / 
Jo 



dk 
~k 



-1 2 



ji{kr] Q ) + 6 / dr]f'(T],k)ji(k(r]o - 77)) 



Vls 



(32) 



where ji(z) is the spherical Bessel function and /(?y, k) is the evolution factor of the gravitational potential defined 
as $(77, k) = $(0, k)f(r), k). The integral over the jfikrjo) term corresponds to the primary SW effect on the last- 
scattering surface, and the integral over the remaining terms corresponds to the ISW effect. The flat CDM model 
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gives /' = and only the primary SW effect, with Kf = 1. The flat ACDM model has no dark energy perturbations 
and gives a scale- independent evolution factor /(ry) which leads to a positive 1SW effect in the low-/ CMB power 
spectrum j^H, rising at the lower-/ (/ = 2,3) end. For the flat QCDM model, the bigger EOS parameter w (> — 1) 
brings out the dark energy earlier and decreases the gravitational potential more, corresponding to a bigger /' value for 
all scales and affecting the magnitude of the CMB power spectrum. The quintessence perturbations grow differently at 
different scales and give a fc-dependent \f'(r], k)\, which is bigger at small scales than at large scales with the transition 
around the horizon size. This means that the lower-/ (/ = 2, 3) multipoles will get suppressed relative to the higher-/ 
ones such that the whole low-/ spectrum will get flattened. These two effects can be seen by numerical calculation 
for the standard ACDM model, constant-^ QCDM model with w — —0.5, and linear- V(Q) QCDM model with 
w(r]o) — —0.5 in Figure [7| This is done with ea. i|32|) approximately for illustrative purpose. Exact results should 
refer to more accurate numerical work like the CMB-FAST code. For non-smooth initial conditions, the perturbation 
evolution can be seen in Table ITU and Table Hvl and corresponds to very unusual cosmological scenarios which we will 
not discuss further. 




2 3 4 5 6 7 

L 



FIG. 7: The low-/ CMB power spectrum due to the Sachs-Wolfe effect in various QCDM models with the smooth initial 
condition. The solid line is for the AC DM model, the dash line for the constant-?/; QCDM model with w = —0.6, and the dot 
line for the linear-V(Q) QCDM model with w(rj ) = -0.5. 



The scale dependence and sensitivity to initial perturbations can also be tested in other ways. The CMB results give 
the primordial $fe spectrum, while the gravitational lensing and other large-scale observations yield $k information 
at low redshift. The ratio of these two, the evolution of gravitational potential, can be directly compared with the 
theoretical plot on different scales shown above for various QCDM models. Furthermore, since the matter density 
perturbation is related to the gravitational potential by eq.JSJ), the evolution of gravitational potential can also be 
inferred from the measurement of the visible matter density spectrum (by SDSS, etc.), and checked with QCDM 
predictions. 

By considering the dark energy perturbations in these QCDM models, we see that they all have a scale dependence 
which is different from the ACMD model, and they have different sensitivity to initial conditions. These properties 
make the dark energy perturbations an effective way to test these models. More precise cosmological observations in 
the future will give better constraints on these QCDM models. 
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